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A new concept of an electromechanical nanodynamometer based on the relative displacement of layers of bilayer 
graphene is proposed. In this nanodynamometer, force acting on one of the graphene layers causes the relative dis- 
placement of this layer and related change of conductance between the layers. Such a force can be determined by mea- 
surements of the tunneling conductance between the layers. Dependences of the interlayer interaction energy and the 
conductance between the graphene layers on their relative position are calculated within the first-principles approach cor- 
rected for van der Waals interactions and the Bardeen method, respectively. The characteristics of the nanodynamometer 
are determined and its possible applications are discussed. 
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I. INTRODUCTION 

Due to the unique electrical and mechanical properties, 
carbon nanostructures (fuUerenes, carbon nanotubes and 
graphene) are considered as promising materials for the use 
in nanoelectromechanical systems (NEMS). Since the con- 
ductance of carbon nanotubes depends on the relative dis- 
placement of nanotube walls at the sub-nanometer scale,i"— 
a set of nanosensors based on such a displacement was pro- 
posed. This set includes a variable nanoresistor,'*'^ a strain 
nanosensor* and a nanothermometer^'** (see Ref. ,9 for a re- 
view). A number of nanotube-based NEMS, such as a nanores- 
onator based on a suspended nanotub&iSiii and a nanoac- 
celerometer based on a telescoping nanotuber^ii^ were sug- 
gested as means for measurements of small forces and accel- 
erations by detection of changes in the system capacitance. 
In addition to zero-dimensional and one-dimensional carbon 
nanostructures, fuUerenes and carbon nanotubes, a novel two- 
dimensional carbon nanostructure, graphene, was discovered 
recently.'^ By analogy with NEMS based on carbon nan- 
otubes, nanodevices based on graphene were proposed.''''^ 
A nanoresonator based on flexural vibrations of suspended 
graphene was implemented.'^ Similar to devices based on the 
dependence of conductance of carbon nanotubes on the rela- 
tive displacement of nanotube walls, NEMS based on the de- 
pendence of conductance of graphene on the relative displace- 
ment of graphene layers can be considered. 

In this paper, we propose a new concept of an electrome- 
chanical nanodynamometer based on the relative displacement 
of layers of bilayer graphene and investigate the operating 
characteristics of this sensor The conceptual design of the 
nanodynamometer is shown in Fig. 1. The operation of the 
nanodynamometer is determined by the balance of only two 
forces: an external force, Fgxt, applied to the movable layer of 
the bilayer graphene, which should be measured, and a force 
of interlayer interaction. Pint- The feedback sensing of the ex- 
ternal force is based on the dependence of the tunneling con- 
ductance G on the displacement of the movable layer under the 
action of the external force. 

The paper is organized in the following way. Calculations of 
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FIG. 1. Conceptual design of the nanodynamometer. The bottom 
graphene layers fixed on the electrodes are indicated as 1 and 2 and 
the movable top graphene layer is indicated as 3. 



the dependence of the force of interlayer interaction on the rel- 
ative displacement of graphene layers and estimations of accu- 
racy of force measurements are presented in Section II. Section 
III is devoted to calculations of the tunneling conductance be- 
tween graphene layers. Our conclusions and discussion of pos- 
sible applications of the nanodynamometer are summarized in 
Section IV. 



II. INTERLAYER INTERACTION OF BILAYER GRAPHENE 

To study the dependence of the force of interlayer interac- 
tion, Fini, on the relative displacement of graphene layers the 
interlayer interaction of bilayer graphene has been investigated 
in the framework of the density functional theory with the dis- 
persion correction (DFT-D). '**''^ The periodic boundary condi- 
tions are applied to a 4.26 A x 2.46 A x 20 A model cell. The 
VASP code^" with the density functional of Perdew, Burke, 
and Ernzerhof^' corrected with the dispersion term (PBE-D)— 
is used. The basis set consists of plane waves with the max- 
imum kinetic energy of 800 eV. The interaction of valence 
electrons with atomic cores is described using the projector 
augmented-wave method (PAW).^^ Integration over the Bril- 
louin zone is performed using the Monkhorst-Pack method^l 
with 24 X 36 X 1 A:-point sampling. In the calculations of the 
potential energy rehefs, one of the graphene layers is rigidly 
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FIG. 2. Structure of a single layer of graphene (see, e.g., Ref. 32) in 
real space (a) and in reciprocal space (b). The elementary unit cell is 
denoted by dotted lines. The hexagon in part (b) is the boundary of 
the first Brillouin zone, ai and ai are the translational vectors, ki and 
k2 are vectors reciprocal to ai and a^; fc,- and ky are projections of k] 
and k2 on coordinate axes. Nonequivalent lattice sites are denoted by 
A and B. 



determined by the condition Fgxt + ^int = 0. This equilibrium 
is stable if the matrix of the second derivatives of the potential 
function U(x,y) is positive definite. Differentiating Eq. ([T), 
we find that the stable equilibrium is possible up to the dis- 
placement l^xil « 0.229acc in the direction corresponding to 
transition from the AB stacking to the SP stacking and up to 
the displacement \6x2\ - acc/4 in the direction corresponding 
to transition from the AB stacking to the A A stacking. The 
maximum forces that can be measured in these directions are 
Fi = 15 pN and F2 = 40 pN per elementary unit cell, respec- 
tively. In Figs. 3 and 4, the regions where the stable equilib- 
rium is possible and is not possible are shown with the solid 
the dashed lines, respectively. 

The upper limit of forces that can be measured using the 
graphene-based nanodynamometer is proportional to the over- 
lap area of the graphene layers and is given by F,^^^ ^ F2NC, 
where A^g is the number of the elementary unit cells in the 



shifted parallel to the other. Account of structure deformation 
induced by the interlayer interaction was shown to be inessen- 
tial for the shape of the potential relief for the interaction be- 
tween graphene-like layers, such as the interwall interaction 
of carbon nanotube s^^'^^ and the intershell interaction of car- 
bon nanoparticles.^^^ The DFT-D calculations show that the 
ground state of bilayer graphene corresponds to the AB stack- 
ing (Bemal structure) with the interlayer spacing 6Z = 3.25 A 
and the interlayer interaction energy -50.6 me V/atom. The in- 
teraction of a single carbon atom in the graphene flake with the 
graphite surface was described using the simple approxima- 
|.Jqj^29,30 containing only the first Fourier components. Based 
on that expression, the interlayer interaction energy U{Sx, 6y) 
as a function of the relative displacements Sx and 6y of the 
layers along the axes x and y chosen along the armchair and 
zigzag directions, respectively, at the equilibrium interlayer 
spacing can be roughly approximated in the form — 



{/, meV 

38.2 



C/= f/i 1.5 + cos 



1 1 1 .5 -H cosl 2kt6x 1 - 

- 2 cosl A:,:^x - - j cos(ky6y) j -1- Uo, (1) 

where ky = 27r/( VSocc), k^ = /ty/ Vs, ace = 1.42 A is the 
bond length of graphene (see Fig. 2), 6x - and ^y = at the 
AB stacking. The parameters L'^o - -101.18 meV and U] = 
8.48 meV (per elementary unit cell) are fitted to reproduce the 
potential energy relief of bilayer graphene. The relative root- 
mean-square deviation 6U/U] of approximation ([TJ from the 
potential energy relief obtained using the DFT-D calculations 
is found to be 6U/U\ = 0.043. The potential energy relief 
calculated using approximation ^ is shown in Fig. 3. 

For simplicity, we restrict the analysis of operation of the 
nanodynamometer by the case where the external force Fext is 
directed along the x (armchair) direction, i.e. along the path 
between adjacent energy minima. The dependence of the in- 
terlayer force on the displacement of the movable layer in this 
direction can be calculated using approximation dU, 



dU_ 
dx 



= lUxkhrnhk^dx - y) - 
- sinl kx6x j co^{ky6y) I. (2) 



This dependence is shown in Fig. 4a. Under the action of 
the external force Fext, the equilibrium position of the layers is 
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FIG. 3. Calculated interlayer interaction energy U of bilayer graphene 
per elementary unit cell as a function of the relative position Sx and 
^V of the layers. The energy is given relative to the global energy 
minimum Uq. SP is a saddle point in the potential relief. The re- 
gions where the stable equilibrium is possible and not possible on the 
displacement of the movable layer along the armchair direction are 
shown with the solid and dashed lines, respectively. 
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FIG. 4. (a) Calculated force Fi„t of the interlayer interaction acting 
on the movable layer of bilayer graphene per elementary unit cell as 
a function of the relative displacement Sx of the movable layer in the 
X (armchair) direction obtained using approximation ([T]. The regions 
where the stable equilibrium is possible and not possible are shown 
with the solid and dashed lines, respectively, (b) Calculated tunneling 
conductance G/Gab between the layers as a function of the relative 
displacement Sx of the movable layer in the x direction. 



overlap area. For example, for the overlap areas of 10^ and 
10"* nm^, the maximum forces that can be measured are 76 nN 
and 7.6 yuN, respectively. The accuracy of the force measure- 
ments is limited by thermal vibrations of the graphene layers. 
The amplitude of these vibrations can be estimated as 



{x^h 



kuT Id^U 



Ng \dx^ 



(3) 



where k^ is the Boltzmann constant, T is temperature and 
d^Uldx^ is the second derivative of the interlayer interaction 
energy with respect to the displacement of the layers along the 
armchair direction at the energy minimum. The latter quan- 
tity is found to be equal d'^UIdx^ = 3Uikl = 55.3 meV/A^ 
per elementary unit cell from Eq. ([T}. The relative error of the 
force measurements can be estimated as the ratio of the am- 
plitude of the thermal vibrations to the maximal displacement 
of graphene layers where the stable equilibrium is possible. 
At liquid helium and room temperatures, these ratios equal 

y/{x^)T/Sx2 = 0.23/ VA^ and yf{x^/6x2 = l.9/^/NE, re- 
spectively. So for the overlap area of 100 nm^, these quantities 
are 0.005 and 0.044, respectively. It is seen that the relative 
error of the force measurements decreases with increasing the 
overlap area of the graphene layers. 



X — (;^(r - Rf ) ± x(r - R^ " d)). (5) 

and the same formula for ^Ptop- Here A^g is the number of the 
elementary unit cells of graphene, d is the vector between two 
non-equivalent carbon atoms (A and B) in the elementary unit 
cell, d - ace', signs + and - correspond to n- (bonding) and 
n*- (antibonding) orbitals in graphene, respectively, R^"' is the 
radius vector of the ^-th unit cell of the bottom graphene layer, 
r is the radius vector The function ;^'(r) is a Slater 2/?^ -orbital 



x(r) = { — \ zexp 



-^-^x^ +y^ +z^\. 



(6) 



where^^ ^ = 1.5679/flB and z is the axis perpendicular to 
the graphene plane; ab = 0.529 A is the Bohr radius and 
r - -\Jx^ +y^ + z} is the magnitude of the radius vector r from 
carbon atom center 

Let us substitute the wave function (|5]l into Eq. (|4]i. The 
product ^i^bot^^top in Eq. ^ can be rewritten as 

= ^ Xi exp(-/kbo,R^°')Or(r - R^°') ±;r(r - Rj°' - d)) x 



III. CONDUCTANCE OF BILAYER GRAPHENE 



Wg 



V 2^ exp(/ktopR7)0r(r - r;?) ±x{r - R"" - d)) 



Let us show that tunneling conductance G between graphene 
layers changes considerably with the relative displacement 
of the layers and therefore measurements of the conductance 
G can be used to determine this displacement. We use the 
Bardeen method,^^ which was previously applied for calcu- 
lation of the tunneling conductance between walls of double- 
walled carbon nanotubes.-^ It is known^"* that the tunneling con- 
ductance is proportional to the sum of squares of the ampli- 
tudes of the tunneling transition (tunneling matrix elements) 
for all electron states at both sides of the tunneling transi- 
tion. This approach was used previously to study the electronic 
structure^^ and conductance^* of twisted two-layer graphene 
system. Here we use such an approach to calculate the rela- 
tive changes of tunneling conductance G between the layers at 
their relative displacement from the ground state correspond- 
ing to AB stacking. 

In the framework of the Bardeen's formalism^^ the ampli- 
tude of the tunneling transition between states of the bottom 
('I'bot) and top (Ttop) layers of bilayer graphene is given by 



/Vg Wg 



11 **^boli''iop 

bot.top 



= 1^ f (^bcV^top - ^topV>p;„,) ofs, (4) 



where S is the overlap area between the graphene layers, kbot 
and ktop are two-dimensional vectors in the reciprocal space 
of the graphene lattice corresponding to the bottom and top 
layers, mo is the electron mass in vacuum, h - hjln is the 
Planck constant. 

In the tight-binding approximation for vectors kbot (or k^p) 
near the comers (/T-points,-'^ K — (In/iSacc), 27r/(3 V3acc)) 
and K' = (27r/(3acc), -27r/(3 V3acc))) of the Brillouin zone 
(Fig. 2b), the wave function of the bottom graphene layer takes 
the form^^ 



= xrp 2j Zj exp(/(k,opR~' - kbotR^^°') + /k^opAR,,) x 

^° S=l li=l 

X Or(r') ±x(r' - d))V(x(r' - AR„) ±x(r' - AR, - d)), (7) 

where r' = r - R^°'. In Eq. (|7]i the coordinates of unit cells of 

the top layer R'°'' - Rp°' + AR/, are expressed via the coordi- 

nates of unit cells of the bottom layer R^°' and displacements 
of unit cells of the top layer AR/,. 

Since all unit cells of graphene are identical, only one cell 
of the bottom layer can be considered in the calculation of 



bot,top 



MJ'^T.J"''- It should also be taken into account that for large 



A^G s> 1 the following relation is satisfied 

, No 

— Y^ exp(/ktopRj) exp(-!kbotRg) = ^k,„„k,„ 



A^G 



(8) 



g=i 



where 5k(,„,.k,„ is the Kronecker symbol. 

Taking into account Eqs. (|7]i and dS), the Eq. (|4]i takes the 
form 

^k, , k,„p ^ y exp(/ktopR^°' - /kbotR^') x 

x ^ exp(!ktopAR/,) I (x(r - AR/,) +x(r - AR/, - d)) x 



h=l 

X VOt(r) ±x(r - d))dS = M,^^,^p<„,k,„^. 



(9) 



Applying 5kb„„k,„p in Eq. Q, we get only vectors k^p = kbot- 
This yields for ktop - kbot = K (or ktop - kbot = K'; see 
Fig. 2b): 
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M. ,, - M. ' - M, ' > - exp(iKAR/,) 

bot,top bot,top bot.top 9i77 1 / i 7 ^ 



2mo fri 2 



xf(x(r- ARh) ± xir - AR, - d))V(x(r) ± x(r - d)) dS ^ 

"a 

^ ^ exp(!KAR;,)(rA-A; + Ta-b;, + Tb-a;, + 7b-b'^), (10) 

where no - [nARl^.^^/( V3aj/2)] is the number of unit cells 
of the top layer located at the distance in graphene plane 
less than ARmax from the considered unit cell of the bottom 
layer, 'Xa(B)-a' (B' ) are the hopping integrals between atom A 
(or B) in the considered unit cell of the bottom layer and 
atom A'l^ (or BJ^) in the h-th unit cell of the top layer We use 

A/?max - 2ai - 2fl2 - 2 V3flcc and nc - [871/ V3] = 14 in the 
calculations for the both layers (Fig. 2a), taking into account 
that the interactions between atoms lying at longer distances 
change the value of the matrix element by less than 0. 1 %. 
The hopping integrals yA(B)-A'{B') in Eq- ( fTOt are given by 

7A{B)-A'^(B'^) =7p(x,y) 

tP- r \ I d d \ 



dz 



dz 



where the index A{B) denotes atom A (or atom E), ;ifbot = 
X{x - XA(E),y - Fa(B), -5Z/2), ^top = x{x - (^a;(bp + 5x), 
y - (Fa'(B') + 5y),SZI2), Xa(b) and Ya(b) are the coordinates 
of atom A (or atom E) in the elementary unit cell of the bottom 
layer {Xa = 0, Xb = occ, Ya^Yb^ 0), ^a;,(b;,) and iA;,(B;,) are 
the coordinates of atoms of the top layer for bilayer graphene 
in the ground state (AB stacking), and dZjl - 1.625 A is 
the half of the interlayer distance, p - {Xa'(b') + 6X - Xa{B), 
Ya'(B') + SY - Ya(B), 0) is the projection on graphene plane of 
the vector connecting two selected atoms in bilayer graphene; 
for jA-A' and jb-b' vector p is the projection of AR/,, for ja-b' 
vector p is the projection of AR/, + d, and for jb-a'. vector p is 
the projection of AR/, - d. Analogously to Refs. 25, and^ the 
hopping integral jp depends on the magnitude p of vector p 
(see Fig. 5). The function jp can be approximated (with accu- 
racy within 3%) by an expression jp - ymax exp(-^(p/flcc)^), 
where ymax = 189 meV and ( = 0.8. 

The calculation according to Eq. (fTTT l for the AA stacking 
of bilayer graphene, in which equivalent atoms of the top and 
bottom layers are located opposite to each other, gives the val- 
ues of the hopping integrals ja'-a - Jb'-b - 189.2 meV. The 
expressions (|6]l, ( fTOl i and (fTTT l allow to calculate the tunnel- 
ing matrix element. For the AB stacking (the Bernal structure, 
6x - Q and 6y - 0), the amplitude of the tunneling transition 
between the states ^Fbot and ^top of the bottom and top layers 



of bilayer graphene is found to be \M. 



K,K 

bot,top' 



136 meV, while 



for the AA stacking {6x = -ace, Sy = 0), W^^^^^] ^ 111 meV 
(see Eq. ^). 

The ratio of the tunneling conductance G to the tunnel- 
ing conductance Gab of bilayer graphene at the ground state 
(6x = 0) equals to the ratio W'^Zo/ I\m'^,1,/s.=o determined 
by Eq. ( fTOl l. The dependence of this ratio on the relative dis- 
placement of the layers along the x (armchair) direction is 
shown in Fig. 4b. It is seen that the tunneling conductance 
between the graphene layers strongly depends on their rela- 
tive position at the sub-nanometer scale, similar to the results 
obtained for double-walled carbon nanotubes.'"^ The conduc- 
tance reaches its maximum for the AA stacking, in which 
atoms of the layers are located at the smallest distances to each 
other. The minimum of the tunneling conductance corresponds 
to the SP stacking. Figure 4 shows that the relative displace- 
ment of the graphene layers in the course of the operation of 
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FIG. 5. The hopping integral yp as a function of the magnitude p of 
vector p (see Eq. l illi l) 



the nanodynamometer can result in changes of the tunneling 
conductance G in the relatively wide range from 0.61 Gab to 
1.73Gab- Thus it is seen that the relative displacement 6x of 
the layers (Fig. 4b) and, consequently, the external force acting 
to the layers (Fig. 4a) can be determined by the measurements 
of the electrical conductance between the layers. 

The model that we use to calculate the tunneling conduc- 
tance adequately describes electron tunneling for relative po- 
sitions of the graphene layers in which their atoms are not 
located exactly opposite to each other. In the case when the 
atoms are located exactly opposite to each other, hybridization 
of their wave functions occurs leading to a significant increase 
of the conductance which is now determined not by tunneling 
between the layers but rather by transitions between energy 
bands of the combined electron system of bilayer graphene. 
Therefore our calculations provide only the lower bound es- 
timate of the relative variation of the tunneling conductance 
upon the relative displacement of the graphene layers. Never- 
theless even this estimate is sufficient to demonstrate the feasi- 
bility of force measurements using the proposed design of the 
nanodynamometer. 

Let us also consider the possibility of a nanodynamometer 
based on the relative rotation of graphene layers. At a relative 
translational displacement of the layers from the ground state 
corresponding to the AB stacking, the interlayer interaction en- 
ergy increases and the tunneling conductance between the lay- 
ers increases or decreases (depending on direction of displace- 
ment) identically for all local areas of the overlap. Contrary to 
that case, at a relative rotation of the layers, the interlayer in- 
teraction energy and the tunneling conductance change differ- 
ently for local areas of the overlap. While the interlayer inter- 
action energy increases for any local area since the AB stack- 
ing corresponds to the global energy minimum, the tunneling 
conductance increases for some local areas and decreases for 
the others. As a result, contributions from different local ar- 
eas to the total tunneling conductance compensate each other 
Therefore changes in the total tunneling conductance at the rel- 
ative rotation of the layers are much smaller than such changes 
at the relative translational displacement. Moreover the force 
required for the relative rotation of the layers from the AB 
stacking to the incommensurate state is an order of magnitude 
greater than the force required for the displacement from the 
AB stacking to the SP stacking.^' Thus the scheme of the nan- 
odynamometer based on the relative rotation of the graphene 
layers is less effective than the proposed scheme based on the 
relative displacement of the layers. For the proposed scheme of 
the nanodynamometer, the relative rotation of the layers should 



be avoided, i.e. only forces that do not produce a significant 
torque should be considered. This is the case when the mea- 
sured force acts 1 ) uniformly on all atoms of the upper layer, 2) 
on adsorbents uniformly distributed on the surface or edges of 
the upper layer in the area between the first and second bottom 
layers 3) on a nanoobject placed near the center of the upper 
layer. 



IV. DISCUSSION AND CONCLUSIONS 

We have proposed the concept of the electromechanical nan- 
odynamometer based on bilayer graphene in which the force 
is determined by measurements of the conductance between 
the layers. In this nanodynamometer, the force acting on one 
of the graphene layers causes the relative displacement of this 
layer and related change of the conductance between the lay- 
ers. The calculations of the potential relief of the interlayer in- 
teraction energy within the dispersion-corrected density func- 
tional theory approach showed that the stable equilibrium of 
bilayer graphene is possible if the measured force acting on 
one of the layers along the armchair direction does not exceed 
40 pN per elementary unit cell. The corresponding displace- 
ment of graphene layers lies within 0.36 A. The calculations 
of the tunneling conductance of bilayer graphene using the 
Bardeen method allowed us to estimate that on the relative dis- 
placement of the layers, the tunneling conductance changes by 
at least a factor of 2, which provides the excellent possibility 
to determine the force by the conductance measurements. The 
relative error of the force measurements is determined by the 
relative thermal vibrations of the layers. This error decreases 
with the increase of the overlap of the layers and with the de- 
crease of temperature. 

Let us discuss possible applications of the considered nano- 
dynamometer A molecule or a nanoobject can be adsorbed on 
the top layer of the nanodynamometer in the region where the 
top layer does not overlap with the bottom layers. The mea- 
surements of the force acting on the molecule or nanoobject 
in the presence of an electric or magnetic field would allow to 
determine their polarizability, electric and magnetic dipole and 
quadrupole moments. 

In the pioneering work of Novoselov et al graphene flakes 
were placed on an insulating substrate and brought into con- 
tact with electrodes''* (to create field-effect transistor). A fur- 
ther considerable progress has been achieved in manipulation 
of individual graphene layers. Individual graphene flakes were 
moved on a graphite surface by the tip of the friction force 
microscope.— The possibilities to cut graphene nanoribbons 
with desirable geometrical parameters'**^ and remove individual 
graphene layers in a controllable way for device patterning'" 
were demonstrated. The tunneling conductance can be mea- 
sured for graphene in the way similar to the experiments for 
multiwall carbon nanotubes. '*-'*'' All these give us a cause for 
optimism that the proposed graphene-based electromechanical 
nanodynamometer will be implemented in the near future. 
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